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ABSTRACT 


NOAA-7  AVHRR  data  from  April  1982  was  used  to  perform  global-scale  analysis 
of  aerosol  particle  characteristics.  Mahony's  (1991)  water  vapor  correction  was  incor¬ 
porated  into  the  AVI  IRR  multichannel  satellite  data  processing  technique  used  by  Frost 
(19S8).  Channel  4  Channel  5  brightness  temperature  difference  was  used  as  an  estimate 
of  water  vapor  content  in  an  air  column.  Greatest  measured  water  vapor  content  was 
at  the  equator,  decreasing  toward  the  poles.  Applying  the  correction  reduced  the  aerosol 
particle  size  index,  resulting  in  an  increase  in  aerosol  optical  depths.  In  addition,  the  high 
particle  size  index  over  the  low  latitudes  in  the  southern  hemisphere  noted  by  Frost,  was 
reduced  significantly  after  applying  the  correction.  A  comparison  was  made  between 
Frost's  values  for  the  variable  two-term  Henyey-Greenstein  phase  function  and  those 
derived  with  the  water  vapor  correction.  The  basic  shape  of  the  curve  is  similar  to  that 
determined  by  Frost. 
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I.  INTRODUCTION 


A.  PURPOSE 

The  availability  of  satellite-measured  radiances  provides  an  effective  and  relatively 
inexpensive  means  for  determining  aerosol  particle  distribution  on  a  global  scale.  Satel¬ 
lites  provide  extensive  coverage  of  the  earth's  surface  with  sufficient  temporal  periodicity 
to  allow  for  the  study  of  large-scale  global  effects  while  at  the  same  time  providing 
small-scale  local  measurements  desirable  for  military  applications.  Data  from 
meteorological  satellites  require  careful  processing  in  order  to  distinguish  the  much  lower 
radiances  associated  with  aerosol  optical  depth  variations  (Durkee  19S4)  from  the 
greater  radiances  associated  with  cloud  albedos,  satellites  were  designed  to  measure. 

Improved  aerosol  optical  depth  retrieval  algorithms  have  facilitated  ongoing  efforts 
to  establish  a  global  aerosol  particle  data  base.  Frost  (19S8)  developed  an  algorithm  to 
derive  aerosol  particle  characteristics  from  the  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  on  NOAA  polar  orbiting  satellites.  Aerosol  particle  size  index 
can  be  derived  by  a  comparison  of  the  sensitivities  of  two  adjacent  spectral  windows  and 
so  detect  aerosol  scattering.  Frost  examined  the  performance  of  a  variable  phase  func¬ 
tion  as  an  improvement  to  the  constant  phase  function  used  previously.  Mahony  (1991) 
examined  the  effects  of  water  vapor  on  aerosol  optical  depth  retrieval.  An  estimate  of 
water  vapor  content  in  a  column  of  air  can  be  derived  by  comparing  the  brightness 
temperature  of  two  AV1IRR  windows  in  the  thermal  infrared  region.  He  concludes  that 
this  effect  can  be  compensated  for,  eliminating  an  error  of  15%  in  the  phase  function. 

The  purpose  of  this  study  is  to  incorporate  Mahony  s  improvement  into  the  algo¬ 
rithm  used  by  Frost,  which  did  not  account  for  water  vapor  appreciably  affecting  re¬ 
trieved  radiances.  The  April  19S2  data  examined  by  Frost  will  be  re-evaluated  to 
compensate  for  water  vapor.  The  variable  two-term  Henyey-Greenstein  phase  function 
recommended  by  Frost  will  be  examined  in  light  of  the  new,  water  vapor  corrected  par¬ 
ticle  size  phase  function  and  compared  to  the  values  he  generated. 

B.  MOTIVATION 

Aerosols  consist  of  a  variety  of  solid  or  liquid  material  dispersed  in  the  atmosphere; 
examples  include  volcanic  ash,  smoke,  and  pollen.  These  particles  vary  in  their  distrib¬ 
ution  and  concentration  in  the  atmosphere,  which  can  impact  local  and  global  radiation 
budgets.  The  extreme  concentrations  of  smoke  particles  resulting  from  the  oil  fires  in 
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Kuwait  created  a  layer  that  reduced  daylight  to  night  in  the  region.  Even  modest  con¬ 
centrations  have  demonstrated  an  impact  in  models  of  the  earth's  radiation  budget 
(Weiss  et  al.  1974;  and  Coakley  1976).  A  Soviet  study  of  longwave  radiation  in  the 
tropical  Atlantic  ocean  (Zaitseva  1976)  linked  short-term  temporal  changes  in  temper¬ 
ature  to  changes  in  aerosol  particles. 

Charlson  et  al.  (19S7)  proposed  that  the  major  source  of  cloud  condensation  nuclei 
(CCN)  over  remote  unpolluted  oceans  appears  to  be  a  waste  product  of  phytoplankton, 
aqueous  dimethylsulphide  gas  (DM S).  A  complex  chain  of  chemical  and  biological 
processes  result  in  DMS  passing  through  the  sea  surface.  Once  in  the  atmosphere,  it  is 
oxidized  to  form  sulfate  aerosol  particles  which  act  as  CCN.  Varitions  in  phytoplankton 
populations  with  the  associated  changes  in  CCN  particles  could  alter  the  earth's  heat 
budget. 

Of  direct  importance  to  the  military  is  the  impact  of  aerosols  on  electro-optical  de¬ 
vices.  The  .4  to  1.0  fim  waveband  is  considered  the  optimum  window  for  lasers  as  it 
avoids  ozone  absorption  at  shorter  wavelengths  and  water  vapor  absorption  at  longer 
wavelengths  (Blocmbergen  et  al.  19S7).  Aerosol  scattering  is  a  significant  factor  in  this 
same  waveband. 
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II.  THEORY 


A.  DIFFUSE  RADIANCE 

Knowing  that  optical  depths  are  small  over  the  ocean  and  that  single  scattering  is 
dominant  (Pfeil  19S6),  and  assuming  an  optically  thin  atmosphere.  Frost  (19SS)  used  the 
following  simplified  version  of  the  Radiation  Transfer  Equation  (RTE): 

2.1 


where, 

L  =  diffuse  radiance, 

4>  =  azimuth  angle, 

oj0  =  single  scatter  albedo, 

F0  =  solar  irradiancc, 

H  =  cos  0,  (0  =  satellite  zenith  angle), 
p  =  scattering  phase  function, 

O  =  single  scatter  angle, 

5  =  optical  depth. 

Griggs  (1975)  and  Durkee  et  al  (1986)  demonstrated  this  linear  dependence  of  radiance 
on  optical  depth. 

There  are  three  primary  sources  of  satellite-detected  radiance;  aerosol  scatter  ( LA ), 
Rayleigh  scatter  ( LR ),  and  sea  surface  reflectance  (Ls).  This  relationship  can  be  expressed 
as 


L  =  LA  +  LR  +  Ls  2.2 

.  Rayleigh  scatter  occurs  when  the  size  of  the  aerosol  particle  is  much  smaller  than  the 
wavelength  of  the  incident  radiation  and,  though  strongly  wavelength  dependent,  it  does 
not  vary  spatially.  Sea  surface  reflectance  is  very  small  for  red  wavelengths,  the  albedo 
is  0.5%,  and  zero  for  radiation  with  wavelengths  longer  than  0.7/<m,  except  where  the 


sun-earth-satellite  geometry'  results  in  sunglint  (Ramsey  1968).  Therefore,  by  avoiding 
data  contaminated  by  sunglint,  total  radiance  can  be  approximated  by  the  sum  of  the 
radiances  due  to  Rayleigh  and  aerosol  scatter  (L  *  LA  +  LR).  By  subtracting  the  compo¬ 
nent  due  to  Rayleigh  scatter,  aerosol  particle  radiance  (LA)  will  remain: 

La  =  L(6,  M,  4>)~Lr~  2.3 

where  SA  is  the  aerosol  particle  optical  depth. 

With  Rayleigh  scatter  forming  a  baseline  radiance,  aerosol  scatter  is  the  dominant 
source  of  upwelled  radiance  variations  (  AL  ~  A LA  ).  Eqn  2.3  indicates  satellite  measured 
radiance  is  directly  related  to  the  changes  in  the  aerosol  particle  characteristics  of  cu0, 
p(Q)  and  5A  .  It  is  apparent  that  a  knowledge  of  p(0)  variations  is  required  to  determine 
a  more  accurate  <5,,  from  LA. 

B.  OPTICAL  DEPTH 

Optical  depth  is  a  measure  of  the  scattering  and  absorption  material  present  in  the 
atmosphere.  It  is  determined  by  vertically  integrating  the  extinction  coefficient  (o,,,), 
from  the  earth's  surface  to  the  top  of  the  atmosphere  (II),  expressed  by 

<5=f  oex,dz.  2.4 

Jo 

Extinction  is  due  to  both  absorption  (aofcJ)  and  scattering  (< ?,„,)  of  radiant  energy;  how¬ 
ever,  in  the  marine  environment,  absorption  can  be  neglected,  leaving  extinction  due  to 
scattering  alone  (<?„,  ^  cr,ccr)  such  that 


rtf 

5  —  I  oscatd^.  2.5 

Jo 

The  amount  of  scattering  is  determined  by  the  size  of  the  particles,  their  scattering 
efficiency,  and  the  size  distribution  of  particles  in  the  column  given  by 


a  scat 


{m,  r)n{r)dr. 


2.6 


where. 
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nr 3  =  particle  cross  sectional  area, 

Qsca,{m,  X,  f)  =  scattering  efTiciency, 
m  =  complex  index  of  refraction, 

X  =  wavelength, 
r  =  particle  radius, 

n(r)  =  number  distribution  of  particles. 

Aerosol  particles  are  most  efficient  in  scattering  radiant  energy  at  wavelengths  near 
their  radios,  becoming  less  efficient  as  the  difference  between  the  two  increases.  Any 
changes  in  the  size  distribution  will  produce  a  corresponding  change  in  SA  and  in  the 
observed  radiance  LA.  The  variations  in  radiance  between  two  wavelengths,  specifically 
red  and  near-infrared,  can  be  quantified  by  the  ratio  (Durkee  et  al.  1991) 

(La),,:,  ~  ' 

This  ratio  is  called  the  particle  size  index  and  is  sensitive  to  variations  in  the  slope 
of  the  particle  size  distribution,  increasing  as  the  slope  increases.  Small  particles  have  the 
greatest  influence  if  the  slope  is  steep,  while  large  particles  are  most  important  for  more 
gradual  slopes.  Generally,  (LA)red  ,  due  to  small  particle  scattering,  is  greater  than  ( LA)n,r 
resulting  in  an  Sn  ratio  larger  than  1  (Mahony  1991). 

C.  WAT  rR  VAPOR  ABSORPTION 

The  dominant  influence  on  retrieved  radiances  in  the  visible  and  near-infrared 
wavelengths  is  aerosol  scattering.  Though  water  vapor  effects  are  not  as  significant  as 
aerosol  scatter,  they  can  result  in  an  over-estimation  of  S,,  which  could  be  interpreted 
as  a  change  in  aerosol  characteristics  or  on  the  other  hand,  could  mask  a  significant 
change.  As  the  desired  result  is  the  estimation  of  aerosol  particle  characteristics  from 
satellite  observations,  inclusion  of  water  vapor  effects  will  improve  the  accuracy  of  re¬ 
sults. 

In  Figure  1,  the  effect  of  water  vapor  can  be  seen  in  the  transmittance  (t)  for 
AVHRR  channels  1  and  2.  Note  that  channel  1  r  is  greater  than  channel  2,  indicating 
less  absorption  is  taking  place.  Therefore,  L2  would  tend  to  decrease  more  due  to  ab¬ 
sorption  by  water  vapor  than  L,  leading  to  the  over  estimation  of  Sl2. 
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Finnic  I.  Atmospheric  transmittance  due  to  water  vapor:  Transmittance  is 
shown  with  NO  A  A  &  AV1IRR  channels  1  and  2  response  functions 
(Mahony  1991). 


1).  BRIGHTNESS  TEMPERATURE  DIFFERENCE 

An  estimate  of  the  water  vapor  content  in  an  air  column  can  be  made  using  the 
AV11RR  10  to  13  /mi  window  region.  In  channels  4  and  5,  at  wavelengths  of  11  and  12’ 
/mi  respectively,  absorption  by  watei  vapor  is  dominant  while  aerosol  scattering  is  rela¬ 
tively  insignificant.  I  he  benefits  of  using  these  channels  is  the  colocation  and  temporal 
linkage  with  channels  1  and  2,  providing  a  water  vapor  snapshot  along  the  same  path 
as  the  aerosol  optical  depth  (5A)  measurements.  Although  not  the  ideal  sensor  windows 
lor  water  vapor  content,  their  use  is  justified  within  the  accuracy  of  the  approximations 
made  (Mahony  1991). 
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The  brightness  temperature  (T)  is  determined  from  the  radiance  received,  due  to 
direct  transfer,  from  the  ocean  surface  to  the  satellite  sensor.  Therefore,  the 
transmittance  (r)  is  directly  proportional  to  the  retrieved  T.  Only  absorption  reduces  the 
radiance  reaching  the  sensor;  scattering  is  relatively  insignificant. 

Dalu  et  al.  (1981)  proposed  a  method  for  determining  sea  surface  temperature  using 
brightness  temperature  difference  corrected  for  air  column  water  vapor  content.  After 
manipulation  of  Dalu's  formula,  the  vertically  integrated  water  vapor  density  («•)  is 
provided  by 


w  =  A(Tt  —  r5)  cos  6. 


2.8 


where  A  is  a  constant.  Using  A  =  19600g°A’~'  nr2,  a  scanning  angle  6  =  0,  and  a  variety 
of  temperature  and  humidity  profiles,  Dalu  (19S6)  applied  this  relationship  to  a  radiative 
transfer  model.  The  resulting  line  graph  (Figure  2)  has  a  correlation  coefficient  of 
R  =  0.99,  with  an  error  of  ±  1.5Ag  m~2  in  density. 

Mahony  determined  that  the  water  vapor  influence  on  the  particle  size  index  was  a 
function  of  w  and  Sn{dry),  the  latter  being  the  magnitude  of  Sn  when  w  =  0.  The 
equation  for  the  curve  that  best  fit  his  data  is 


Sn{drv) 


S 12 

1  +  O.U332n/u-- 


2.9 


where  Sl2  and  iv  are  determined  from  measured  radiances. 


E.  VARIABLE  PHASE  FUNCTION 

Based  on  Frost's  work,  a  variable  single  scatter  phase  function  was  used,  having 
been  determined  to  be  more  accurate  in  its  optical  depth  estimation.  Higher  ratios,  which 
are  due  to  smaller  particles,  have  more  backscatter  than  the  lower  ratios,  due  to  larger 
particles.  This  provides  the  basic  shape  necessary  for  the  phase  function  needed  to  esti¬ 
mate  optical  depth. 

The  two-term  Henyev-Greenstein  scattering  phase  function  (Lenoble  1985)  was 
used. 


7 


figure  2.  Total  water  vapor  content  as  a  function  of  brightness  temperature  differ¬ 
ence  (i)alu  1986). 


P  (©)  = 


_ H'd-gf) 

(1  +  S2\  ~  cos  0) 


,  (1  -  ff"Kl -^f) 

(l  +  £?  +  ftcos0)  ’ 


2.10 


where  IT  is  a  weighting  factor,  g,  and  g2  arc  asymmetry  curve  factors,  and  0  is  the  scatter 
angle.  All  three  of  these  factors  arc  functions  of  Sn  written  as 

H'«  1.2 -(0.25  x  Su),  2.11 

ft  =  1.04  -(0.08333  x  S12),  2. 12* 

ft -  1.2- (0.58333  x  S12).  2.13 


from  the  equations  2.11  through  2.13,  the  phase  function's  dependence  on  Su  is  appar¬ 
ent.  The  marine-type  phase  function  describes  large  particle  distribution,  low  Sl2  ratios; 
while  the  rural-type  phase  function  follows  small  particle  distribution,  high  St!  ratios. 
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As  the  ratios  change,  the  factors  change,  resulting  in  a  varying  phase  fuction  from  one 
similar  to  a  marine  phase  function  to  one  similar  to  a  rural  phase  function.  Figure  3 
shows  some  of  the  variable  phase  functions  />(©)  calculated  by  Frost  (19S8)  from  the 
indicated  Sl2  values. 
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Figure  3. 


Aerosol  particle  size  index  and  associated  phase  functions:  (Frost  I9SS). 
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III.  PROCEDURES 


A.  DATA  AND  PROCESSING 

The  data  examined  was  collected  by  the  NOAA  7  polar  orbiting  satellite  whose  av¬ 
erage  altitude  was  833  km.  This  sun-synchronous  platform  crossed  the  equator  at  ap¬ 
proximately  0230  and  1430  local  standard  time  (north-  and  south-bound,  respectively). 
Of  the  five  spectral  bands  or  channels  on  the  AVHRR  instrument;  channels  1,  2,  4,  and 
5  were  used  in  measuring  upwelling  radiance.  Table  1  provides  the  bandwidths  for  all  the 
channels  (Jensen  1986)  plus  a  weighted  average  of  the  solar  irradiance  (fj,)  over  the 
bandwidth  used  for  channels  1  and  2  calculations. 


Table  1.  NOAA  7  AVHRR  CHANNEL  BANDWIDTHS. 


Channel 

Bandwidth 

Solar  Irradiance 

1 

0.58-0.68  uni 

154.3536  ntlVenr1 

2 

0.725- 1 . 10  tint 

103.09US  /nil  car 2 

3 

3.55-3.93  nm 

4 

10.5-11.3  uni 

5 

1 1.5-12.5  urn 

The  archives  of  the  National  Center  for  Atmospheric  Research  (NCAR)  in  Boulder, 
Colorado,  provided  AVHRR  Global  Area  Coverage  (GAC)  data  for  two  periods  in  April 
1982;  the  5th  -  10th  and  the  20th  -  25th.  April  1982  data  was  analyzed  for  comparison 
with  the  work  done  earlier  by  Frost  (19S8). 

All  the  data  was  processed  on  the  Interactive  Environmental  Digital  Analysis  Lab¬ 
oratory  (IDEA  LAB)  located  at  the  Naval  Postgraduate  School,  Monterey,  CA.  The 
IDEA  LAB  produced  enhanced  images  and  graphs  from  the  processed  data. 
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B.  ANALYSIS 


The  data  analysis  scheme  was  based  on  that  used  by  Pfiel  (1986)  and  Frost  (1988). 
The  water  vapor  algorithm  used  by  Dalu  (1986)  was  incorporated  to  account  for  ab¬ 
sorption  by  water  vapor,  as  recommended  by  Mahony  (1991).  Also  accounted  for  in  the 
algorithm  were  changes  to  the  solar  irradiance  due  to  variations  in  the  sun-earth  dis¬ 
tance.  During  preprocessing,  pixels  located  at  latitudes  greater  than  70°  or  those  over 
land  were  discarded. 

To  remove  sunglint,  an  empirical  method  using  viewing  geometry  was  applied.  This 
method,  a  weighted  sum  of  relative  zenith  and  azimuth  angles  for  the  sun  and  satellite, 
identified  approximately  the  same  contaminated  areas  as  the  Avian  sunglint  routine, 
which  is  based  on  the  Cox  and  Munk  (1954)  method. 

Ozone  absorption  effects  occur  in  the  red  and  near  infrared  wavelengths  with  the 
greatest  absorption  at  a  wavelength  of0.69^tfj  (Fleagle  and  Businger,  1980),  which  is  in 
the  channel  1  window.  This  effect  was  inadvertently  compensated  for  twice  by  Frost.  The 
error  has  been  corrected  for  this  study. 

Each  pixel  was  then  examined  for  clouds  as  only  clear  pixels  were  to  be  processed. 

Cloud  removal  was  performed  by  the  same  method  used  by  Frost: 

1.  Test  for  bad  pixels  —  If  channel  2  or  channel  4  radiance  values  are  less  than  zero, 
after  Rayleigh  radiance  was  removed,  then  the  pixel  is  classified  as  "bad"  and  not 
processed. 

2.  Test  for  high  clouds  --  A  pixel  was  considered  to  contain  high  clouds  if  its  channel 
4  brightness  temperature  was  less  than  213° K. 

3.  Test  for  low  clouds  —  A  pixel  was  considered  to  contain  low  clouds  if  its  channel 
2  albedo  was  greater  than  40%. 

4.  Preliminary  test  for  clear  pixels  —  A  channel  1  channel  2  albedo  was  calculated  and 
if  that  ratio  was  greater  than  1.5,  the  pixel  was  tentatively  classified  as  clear. 

5.  Final  test  for  clear  pixel  —  The  channel  2  albedo  of  the  pixel  tentatively  identified 
as  clear  was  compared  with  the  four  pixels  above,  below,  to  the  right  and  to  the 
left.  All  the  pixels  must  fall  with  in  five  counts  of  each  other  for  the  pixel  to  be 
classified  as  clear,  otherwise  it  was  counted  as  partly  cloudy.  A  count  is  equivalent 
to  0.1069%  albedo  (Kidwell  1986). 
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The  clear  pixels  were  then  processed  for  a  variety  of  parameters.  The  list  is  basically 

the  same  as  Frost's,  the  most  significant  difference  is  the  inclusion  of  two  parameters 

with  the  water  vapor  correction,  the  terms  include  the  word  "dry"  in  their  name. 
Channel  1,  Channel  2  ratio, 

Dry  Channel  1  Channel  2  ratio, 

Aerosol  optical  depth, 

Dry  Aerosol  optical  depth, 

Average  phase  function  of  Channels  1  and  2, 

Channel  4,  Channel  5  brightness  temperature  differences, 

Number  of  clear  pixels, 

Number  of  partially  cloudy  pixels. 

P  parameter. 

The  P  parameter  is  used  to  examine  the  phase  function  calulated  by  Frost. 

The  phase  function  p{9,  S,3)  is  an  important  factor  in  the  determination  of  aerosol 
optical  depth.  Any  changes  made  in  the  phase  function  algorithm,  such  as  the  addition 
of  a  water  vapor  correction  to  Sa,  will  produce  corresponding  changes  to  the  calculated 
optical  depths.  A  comparison  of  this  modified  phase  function  to  that  determined  by 
Frost  should  provide  confirmation  that  the  changes  made  produce  physically  sound  re¬ 
sults  and  add  to  the  knowledge  of  p{6,  Su). 

A  new  parameter  was  defined  using  the  simplified  RTF  (2.1)  that  is  calculated  from 
radiance  data  and  can  be  used  to  derive  p{9)  and  SA: 

4pL, 

p  =  vP(e,  sn)dA.  3.i 

/o 

Using  the  fact  that  integrating  the  phase  function  over  a  sphere  is  unity,  the  following 
form  can  be  used  to  determine  p{9,  S12): 

Jp(0,S12)</Q  =  ^Ap{6,Sn)6AdQ.  =  jAPdQ  =1,  3.2 

where  A  is  equal  to  the  inverse  of  the  optical  depth.  It  is  this  equation  from  which  the 
phase  function  will  be  calculated  for  comparison  to  Frost's  curves. 
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All  of  the  results  were  stored  and  averaged  in  1°  by  1°  boxes;  each  data  point  is  one 
box.  The  number  of  pixels  varied  from  box  to  box,  depending  on  the  number  of  partially 
cloudy  to  cloudy  pixels  found  in  each  square. 


14 


IV.  RESULTS 


A.  INTRODUCTION 

Ongoing  efforts  to  determine  aerosol  particle  characteristics  have  resulted  in  a 
steady  improvement  in  the  algorithm.  Water  vapor  effects  were  shown  to  be  an  impor¬ 
tant  factor  in  aerosol  optical  depth  retrieval  but  have  not  been  previously  applied  to 
satellite  data.  In  this  study  the  latitudinal  variation  in  water  vapor  content  is  initially 
examined  to  give  an  overview  of  where  the  effects  could  be  the  most  significant.  Particle 
distribution  characteristics  arc  examined  next,  both  with  and  without  the  water  vapor 
correction.  The  data  having  the  correction  will  be  refered  to  as  "dry".  Both  the  corrected 
and  uncorrected  particle  size  index  and  aerosol  optical  depth  will  be  examined  and 
compared  over  the  Pacific  Ocean  followed  by  the  latitude  profiles  of  the  averaged  values. 
Finally,  the  shape  of  the  P  parameter  curves  are  compared  to  Frost's  phase  function 
curves. 

B.  COLUMN  WATER  VAPOR 

Figures  4  and  5  show  the  average  water  vapor  content  in  an  air  column  to  provide 
a  general  overview  of  water  vapor  effects  in  the  atmosphere  for  the  periods  of  April  5  - 
10  and  20  -  25,  respectively.  The  genera!  structure  is  that  of  a  maximum  in  the  vicinity 
of  the  equator  with  values  then  decreasing  toward  the  poles. The  drop  is  more  rapid  in 
the  northern  hemisphere  than  the  southern  hemisphere,  the  latter  is  transitioning  from 
its  summer  season  when  temperatures  are  higher  and  more  water  vapor  is  present.  The 
greater  year-round  air  and  sea  surface  temperatures  at  the  equator  allow  more  water 
vapor  to  exist  in  a  column  of  air.  From  this  it  can  be  seen  that  the  greatest  potential  for 
error  due  to  water  vapor  is  in  the  equatorial  region.  An  examination  of  the  particle  size 
index  before  and  after  applying  the  water  vapor  correction  will  be  made  next. 
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Figure  4.  Brightness  temperature  difference  fur  April  5  -  10. 
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C.  PARTICLE  SIZE  INDEX 


Figures  6  and  7  depict  the  uncorrected  S12  values  across  the  Pacific  Ocean  for  April 
5-10  and  20  -  25,  respectively.  A  broad  region  of  high  values  extends  across  the  low 
latitudes  in  the  southern  hemisphere  in  both  figures,  indicating  the  presense  of  a  greater 
number  of  smaller  particles.  Low  values  indicating  larger  particles  can  be  seen  in  Figure 
6  across  the  northern  hemisphere.  In  Figure  7  there  are  two  noticable  maxima;  one  in 
the  vicinity  of  the  New  Hebrides  islands,  the  second  west  of  South  America.  The  north 
Atlantic  appears  to  be  a  region  of  small  particles  except  along  the  east  coast  of  the 
United  States. 

In  Figures  6  and  7,  bands  cross  each  figure  from  north-northwest  to  south- 
southeast.  These  are  a  result  of  the  moving  of  the  sensor's  aperature  to  the  left  and  right 
of  the  satellite's  path  to  view  the  side  regions.  Aerosol  scatter  is  anisotropic,  favoring 
forward  and  backward  scatter.  A  satellite  is  looking  at  backscatter  when  it  looks  directly 
at  the  earth.  But  as  the  sensor  looks  to  the  side,  it  is  no  i^ger  viewing  backscatter,  but 
the  smaller  scatter  associated  with  smaller  scattering  angles.  This  results  in  a  decrease 
in  the  phase  function  and  a  corresponding  increase  in  the  optical  depth.  Future  work 
will  need  to  correct  for  this  effect. 

Figures  8  and  9  show  a  comparison  of  the  particle  size  index  profiles  for  Sn  and 
Sn(dry)  for  early  and  late  April,  respectively.  The  Su  profiles  in  both  figures  show  a 
broad  region  of  high  values,  approximately  2.35,  from  the  equator  to  40°S;  the  values 
then  decrease  toward  the  poles  more  rapidly  in  the  southern  hemisphere  than  in  the 
northern.  The  greatest  drop  can  be  seen  in  the  southern  hemisphere  in  Figure  8  with  a 
minimum  value  of  2.0.  The  minimum  in  Figure  9  is  2.10.  Frost's  values  were  2.0  at  the 
equator  decreasing  to  1.65  moving  southward  toward  the  pole;  a  second  smaller  decrease 
occurred  toward  the  north  pole,  dropping  to  1.9.  His  underestimation  was  due  to  a 
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LATITUDE  PROFILE 


Figure  8.  Comparison  of  particle  si/e  index  profiles  for  April  5  -  It):  the  solid  line 
is  S(i,  and  the  dashed  line  is 
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Figure  9 
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Comparison  of  particle  size  index  profiles  for  April  20  -  25: 
is  5,j,  and  the  dashed  line  is  Su(Jry). 
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In  early  April,  Figure  10,  distinct  high  value  areas  remain  in  the  Gulf  of  Carpentaria 
and  along  the  coast  of  China.  Late  April  (Figure  11)  shows  two  areas  of  higher  values 
in  the  northern  Atlantic  Ocean  and  along  the  west  coast  of  North  America.  The  maxi¬ 
mum  west  of  South  America  in  Figure  7  is  easier  to  distinguish  in  Figure  1 1  and  has 
lower  values.  The  maxima  off  the  Americas  appear  to  be  associated  with  upwelling  re¬ 
gions  that  are  biologically  very'  productive,  which  could  be  a  source  for  the  greater 
number  of  small  aerosol  particles. 

In  Figures  8  and  9,  it  is  easy  to  see  the  decrease  in  size  index  values  between  S, 3  and 
Sl2(dry).  At  the  equator  Sn(dry)  is  0.5  less  than  S12  in  early  April,  0.4  less  later  in  the 
month.  From  the  equator  to  the  poles,  the  difference  between  S,2  and  S^{dry)  decreases 
to  approximately  half  the  values  at  the  equator  in  both  figures.  In  both  instances  the 
overall  shape  of  the  SJdry)  curve  shows  a  minimum  near  the  equator  gradually  in¬ 
creasing  as  you  move  toward  the  poles,  which  is  the  reverse  of  Frost's  results  where 
values  decreased  toward  the  poles.  Though  a  more  uniform  profile  was  expected,  the 
peaks  that  exist  in  the  corrected  profiles  of  both  figures  correspond  to  the  latitudes  of 
the  maxima  depicted  in  Figures  10  and  1 1. 

E.  OPTICAL  DEPTH 

The  views  of  the  Pacific  in  Figures  12  and  13,  show  the  southern  hemisphere  with 
low  optical  depths  and  the  northern  hemisphere  with  very'  high  optical  depths.  The 
largest  feature  is  the  wide  area  of  high  values.  >.45  to  0.50,  that  occurs  along  the 
Mexican.  Central  American  west  coast  extending  westward  into  the  Pacific.  This  feature 
is  the  result  of  the  El  Chicon  volcano  (17.3°A\  95.2°IF)  erupting  on  a  number  of  occa¬ 
sions  during  the  period  of  28  March  to  6  April.  The  high  values  in  the  Gulf  of  Mexico 
are  also  the  result  of  El  Chicon.  The  movement  to  the  east  is  possibly  due  to 
stratospheric  winds  carrying  particles  ejected  to  those  heights  (Frost  1988).  In  early 
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figure  14.  Comparison  of  optical  depth  values  for  April  5  -  It):  solid  lines  are 
values,  dashed  are  dA(ilry)  values. 
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Figure  15.  Comparison  of  optical  dcpfli  values  for  April  20  -  25:  solid  lines  arc  54 
values,  dashed  are  «5 A(iiry)  values. 
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changes  occur  in  areas  where  the  number  of  aerosol  particles  are  sufficient  to  affect  the 
optical  depth. 

G.  PHASE  FUNCTION 

Figure  18  shows  the  P  parameter  as  a  function  of  scatter  angle  for  three  different 
particle  size  indices.  A  comparison  to  Frost's  plot  of  the  Henyev-Greenstein  phase 
function  (figure  3)  indicates  that  the  basic  shape  of  the  P  parameter  curves  are  similar 
to  that  shown  by  Frost.  A  smaller  particle  size  index  has  a  larger  value  in  the  forward 
direction  (between  0°  and  60°  )  than  larger  index  values.  The  larger  the  particle  size  index 
the  closer  the  curve  approaches  that  of  Rayleigh  scatter,  which  is  expected  as  smaller 
particles  dominate.  Where  the  scatter  angle  is  between  90°  to  140°,  the  P  parameter  ap¬ 
proaches  a  minimum  value  for  each  curve.  As  the  angle  increases  to  180°  the  curves 
show  an  increase  in  P  values. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Frost  showed  that  his  technique  for  processing  AVHRR  data  was  effective  as  a  tool 
in  the  analysis  of  aerosol  particle  characteristics.  Though  significant  improvements  were 
made  with  his  study,  sources  of  error  have  since  been  identified  requiring  futher  im¬ 
provements  to  his  algorithm.  The  determination  that  water  vapor  was  a  significant 
source  of  error  suggested  a  re-evaluation  of  the  April  19S2  data  in  order  to  assess  water 
vapor  impact  on  aerosol  optical  depth  retrieval. 

The  results  of  this  study  have  shown  that  inclusion  of  a  water  vapor  correction 
produces  a  change  of  approximately  15%  in  the  aerosol  particle  size  index.  Without  this 
correction,  large  regions  of  high  particle  size  index  values,  which  correspond  to  a 
predominance  of  relatively  smaller  particles,  were  actually  regions  of  lower  values 
masked  by  water  vapor  effects.  This  made  it  difficult  to  see  true  areas  of  small  particles 
as  they  were  indistinguishable  from  the  water  vapor  variations. 

Since  the  phase  function  algorithm  is  a  function  of  S]2 ,  the  modification  made  for 
water  vapor  will  affect  the  phase  function.  It  has  been  shown  that  the  shape  of  the 
modified  phase  function  is  similar  to  the  one  determined  by  Frost.  The  improved 
parameterization  cannot  be  directly  compared  to  Frost's  until  the  integration  (3.2)  is 
accomplished. 

The  aerosol  optical  depth  has  also  been  shown  to  be  modified  by  the  water  vapor 
correction  producing  an  approximate  change  of  5%.  The  observed  extent  and  magni¬ 
tude  of  areas  of  high  dry  optical  depths  was  greater  than  the  optical  depths  without  the 
correction. 

Additional  improvements  to  the  algorithm  used  are  still  necessary.  Multiple  scatter 
has  not  been  fully  accounted  for.  The  assumption  used  to  allow  for  single  scatter  only 
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is  valid  for  optical  depths  less  than  about  0.2.  This  value  was  exceeded  in  the  regions  of 
high  optical  depth  associated  with  major  aerosol  particle  producing  events  such  as  the 
El  Chicon  eruptions.  Errors  in  the  values  due  to  multiple  scatter  exist  and  must  be  cor¬ 
rected. 

A  second  area  for  improvement  would  be  in  the  determination  of  Sl2.  Current  pro¬ 
cedures  use  a  ratio  of  channel  1 /channel  2  radiances  to  calculate  the  particle  size  index. 
The  preferred  method  is  to  use  a  ratio  of  channel  1  'channel  2  optical  depths,  which  re¬ 
quires  a  better  knowledge  of  the  phase  function  for  both  channels.  It  is  the  optical  depth 
which  contains  information  on  the  aerosol  particles  and  would  provide  the  best  estimate 
for  Su.  In  this  study  the  phase  function  and  optical  depths  are  at  the  channel  1  wave¬ 
length.  This  improvement  should  eliminate  the  banding  that  occurs  in  the  Pacific  basin 
photographs  which  possibly  hides  features  of  importance. 

Finally,  additional  studies  of  global  aerosol  particle  characteristics  are  needed.  This 
study  was  limited  to  a  single  month.  Studies  exploring  the  changes  that  occur  during  a 
year  and  those  consisting  of  a  number  of  years  by  month  or  season  would  increase  our 
understanding  of  the  role  of  aerosol  particles  in  the  global  climate. 


35 


REFERENCES 


Bloembergen,  N.,  C.  K.  Patel,  and  G.  Pake  (Chairpersons),  1987:  Science  and  Technol¬ 
ogy  of  Directed  Energy  Weapons.  Report  of  the  American  Physical  Society  Study 
Group,  Reviews  of  Modern  Physics,  59,  S1-S202. 

Charlson,  R.  J..  J.  E.  Lovelock,  M.  O.  Andrea,  and  S.  G.  Warren,  19S7:  Oceanic 
phytoplankton,  atmospheric  sulfur,  cloud  albedo  and  climate.  Nature,  325,  655-661. 

Cox,  C.  C.,  and  W.  Munk,  1954:  Measurement  of  the  roughness  of  the  sea  surface  from 
photoeraphs  of  the  sun's  glitter.  Journal  of  the  Optical  Societv  of  America,  44, 
838-850. 

Coaklcy,  J.  A..  1976:  Aerosols  and  the  Earth's  radiation  budget.  Radiation  in  the  At¬ 
mosphere.  H.  J.  Bollc,  ed.,  Science  Press,  Princeton,  NJ,  472-474. 

Dalu,  G.,  C.  Prabhakara,  and  R.  C.  Lo,  1981:  Improved  accuracy  of  the  remote  sensing 
of  sea  surface  temperature.  In  Oceanography  from  Space,  J.  E.  R.  Gower,  ed.,  Ple¬ 
num,  New  York,  109pp. 

Dalu.  G.,  19S6:  Satellite  remote  sensing  of  atmospheric  water  vapor.  International 
Journal  of  Remote  Sensing,  7,  1089-1097. 

Durkee,  P.  A.,  1984:  The  relationship  between  marine  aerosol  particles  and  satellite- 
detected  radiance.  Ph.D.  Thesis,  Colorado  State  University,  Fort  Collins,  CO,  US 
ISSN  0067-0340,  124pp. 

Durkee,  P.  A.,  D.  R.  Jensen.  E.  E.  Hindman,  and  T.  II.  Yonder  Haar,  19S6:  The  re¬ 
lationship  between  marine  aerosols  and  satellite  detected  radiance.  Journal  of 
Geophysical  Research,  91,  4063-4072. 

Durkee,  P.  A.,  F.  Pfeil,  E.  Frost,  and  R.  Shcma,  1991:  Global  analysis  of  aerosol  particle 
characteristics.  Atmospheric  Environment,  25A-11,  2457-2471. 

Fleaglc.  R.  W.  and  J.  A.  Busingcr,  19S0:  An  Introduction  lo  Atmospheric  Physics,  Aca¬ 
demic  Press,  New  York,  NY,  432pp. 

Frost,  E.  M.,  1988:  Global  scale  estimates  of  aerosol  particle  characteristics.  Master's 
Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  1-41. 

Griges,  M.,1975:  Measurement  of  atmospheric  aerosol  optical  thickness  o.er  water  us¬ 
ing  ERTS-1  data.  Journal  of  Atmospheric  Pollution  Control  Association,  25,  622-626. 

Jensen,  J.  R.,1986:  Introductory  Digital  Image  Processing;  A  Remote  Sensing  Approach. 
Prentice-Hall,  Englewood  Cliffs,  NJ,  379pp. 

Kidwell,  K.  B.,  19S6:  NO  A  A  Polar  Orbiter  Data  Users  Guide.  NOAA  NESDIS  National 
Climatic  Data  Center,  Satellite  Data  Service  Division, 145pp. 


36 


Lenoble,  J.,  1985:  Radiative  Transfer  in  Scattering  and  Absorbing  Atmosphere:  Standard 
Computionai  Procedures.,  A.  Deepak  Pub.,  Hampton,  VA,  300pp. 

Mahony,  T.  P.,  1991:  Water  vapor  influence  on  satellite-measured  aerosol  character¬ 
istics.  Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  1-38. 

Pfeil,  F.  R.,  1986:  Developing  a  physical  basis  for  an  aerosol  climatology  of  the  pacific 
ocean.  Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  76pp. 

Ramsev,  R.  C.,  1968:  Studv  of  the  Remote  Measurement  of  Ocean  Color.  Final  Report, 
TRW,  NAS W- 1658.  ' 

Weiss,  R.,  R.J.  Charlson,  A.  P.  Waggoner,  M.  B.  Baker,  D.  Cover,  D.  Thorsell  and  S. 
Yuen,  1976:  Application  of  directly  measured  aerosol  radiative  properties  to  climate 
models.  Radiation  in  the  Atmosphere  H.  J.  Bolle,  ed.,  Science  Press,  Princeton,  NJ, 
469-471. 

Zaitseva,  N.  A.,  1976:  The  spatial-temporal  variability  in  long-wave  radiation  fields  at 
the  GATE-74  grounds.  Radiation  in  the  Atmosphere  .  H.  J.  Bole,  ed.,  Science  Press, 
Princeton,  NJ,  530-532. 


37 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5002 

3.  Chairman  (Code  MR  Hy)  1 

Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey.  CA  93943-5000 

4.  Professor  Phillip  A.  Durk.ee  (Code  MR.'De)  2 

Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

5.  Professor  Patricia  M.  Pauley  (code  MR' Pa)  1 

Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

6.  LT  Mary  B.  Clifford  1 

Fleet  Numerical  Oceanography  Center 

Monterey,  CA  93943-5005 

7.  Commander  1 


Naval  Oceanography  Command 
Stcnnis  Space  Center 
MS  39529-5000 

8.  Commanding  Officer  1 

Fleet  Numerical  Oceanography  Center 

Monterey,  CA  93943-5005 

9.  Commanding  OfTicer  1 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory 

Stennis  Space  Center 
MS  39529-5004 

10.  Director  1 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory 

Monterey,  CA  93943-5006 


38 


1 1 .  Chief  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 

12.  Office  of  Naval  Research 

Naval  Ocean  Research  and  Development  Activity 
800  N.  Quincy  Street 
Arlington,  VA  22217 


39 


